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Abstract 

Predictions for energy levels, production and decay rate of gluinonia, non- 
relativistic boundstates of gluinos, are presented. The potential between color-octet 
constituents is derived in next-to-next-to leading order and one-loop QCD correc- 
tions are derived for the production cross section and the decay rate into gluon jets. 
In addition we evaluate the decay rate into top quarks and into two photons. The 
signal-to-background ratio is estimated for the dominant decay mode and found 
to be around 0.5%. For relatively light gluinos the bound states thus might be 
detectable. 



PACS numbers: 12.38.Bx 12.60.Jv 14.80.Ly 



1 Introduction 



With the turning-on of the LHC in the near future, one of the most important tasks of 
both general-purpose detectors ATLAS and CMS will be the search for physics beyond 
the Standard Model, with one of the most promising candidates being Supersymmetry. In 
the course of this program, which will obviously be adjusted and focused with increasing 
luminosity and with the first indications of "New Physics", the determination of the 
properties of newly detected particles will be a major task. Indeed, only after the quantum 
numbers, the interaction and the mass of a new particle have been measured with high 
precision, this requirement can be considered fulfilled. 

In the present paper we want to concentrate on the gluino, the super symmetric part- 
ner of the gluon, and discuss the possibility to determine its quantum numbers, i.e. its 
spin, color, Majorana nature and its mass, through the investigation of non-relativistic 
boundstates of two gluinos, generally denoted "gluinonia" . The classification of gluinonia 
according to their quantum numbers and their qualitative properties, i.e. spectra and de- 
cay modes, as well as estimates of their production cross section at hadron colliders have 
been studied already more than two decades ago and presented in Refs. [1-3] concentrating 
at that time on gluino masses which are by now experimentally excluded. More recently 
gluinonia have been discussed in Refs. [4,5]. Hadronic transitions between different gluino- 
nium levels and some rare annihilation decays were evaluated in [4]. The production of 
vector and pseudoscalar states at the Tevatron and the LHC was investigated in [5] where 
it was argued that, given sufficiently good yet realistic jet mass resolution, pseudoscalar 
states could be detected for bound state masses up to 3 TeV. A phenomenological study 
of the question, to which extent a signal for color-octet states could indeed be detected ex- 
perimentally, has been performed in [6]. In all these cases, however, the production cross 
section has been calculated in lowest order approximation only and a qualitative phe- 
nomenological potential for the gluino-gluino interaction has been employed. The topic 
has gained renewed interest in connection with Split SUSY [7,8], (see, e.g., Refs. [9,10]). 

The detection of these bound states would provide important information about gluino 
properties, difficult to obtain through other means. Since gluinonia only exist if gluinos do 
not decay too rapidly, say with decay rates less than a few GeV, their observation would 
immediately provide at least an upper limit on T^. The energy spectrum of the gluinonium 
system, e.g. the mass difference between ground state and first radial excitation (or open 
gluino production) would be sensitive to the potential and thus to the color-octet nature of 
the constituents. The same is true for the production cross section which is proportional 
to the square of the wave function at the origin. Gluinos are expected to decay through 
cascades into several jets plus the lightest supersymmetric particle (LSP) which escapes 
detection. For this reason the precise determination of its mass in direct gluino decays is 
limited by the missing mass or energy resolution. 

In contrast, once the mass of the boundstate and the dynamics of the interaction are known 
with sufficient precision, the mass of the constituents is fixed with the corresponding 
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accuracy and in any case subject to a different systematic error. 

The energy levels of gluinonia were estimated already in the early papers employing phe- 
nomenological color-octet potentials, obtained from phenomenological quarkonium poten- 
tials by multiplication with the ratio of the corresponding Casimir operators Ca/Cf = 9/4 
and the gluino mass was interpreted as an effective mass closely related to the pole mass. 
During recent years more refined studies have been performed to describe non-relativistic 
boundstates of color-triplet constituents. This includes the evaluation of the static two- 
loop potential [11-13], first steps toward the evaluation of the three-loop potential [14] and 
the calculation of mass-suppressed terms proportional a 2 s /m [15-19] and a s /m 2 [17,19-21]. 
Mass definitions better suited to control the long distance part of the potential have been 
introduced [22,23], and the evaluation of the spectra through analytic methods has been 
possible, using a perturbative approach to include higher order contributions from the 
potential [24-26]. Last but not least, finite- width effects were incorporated for top-quark 
production at electron-positron [27-33] and hadron colliders [34-36]. These new results 
can be applied also in studies of gluino bound states. 

In order to arrive at reliable predictions for the signals of gluinonium production and 
decay, the inclusion of a realistic QCD potential and next-to-leading order (NLO) correc- 
tions to boundstate production at the LHC is mandatory. Considering and building on 
the significant progress in the perturbative treatment of non-relativistic boundstates, in 
particular quarkonium and positronium, predictions of similar quality can be made for 
boundstates with color-octet constituents. This is the main content of the present paper. 

We start with a brief overview and a discussion of the quantum numbers and qualitative 
properties of boundstates with Majorana constituents and evaluate the QCD potential 
for a color-octet boundstate (Section [2]). Subsequently, in Section [31 we calculate the 
spectrum and the wave function at the origin for the lowest-lying levels. 

These quantities are the necessary ingredients for a realistic prediction of gluinonium 
production and decay rates presented in Section |U Assuming that squarks are rela- 
tively heavy, the dominant decay channel of the pseudoscalar singlet boundstate proceeds 
through two gluons, and it will be a difficult task to distinguish this signal from the 
irreducible background from the two-gluon jet continuum. 

In Section [5] we consider as an alternative the decays into top quarks and two photons, 
which can proceed through virtual quark (q) and squark (q) intermediate states. One 
might, in principle, hope that the extremely clean 77 signal sticks out of the irreducible 
background similarly to the Higgs decay into two photons which has been identified as a 
promising signal significantly superior to the bb or two-gluon mode. Similar considerations 
apply to the tt mode. 

In order that these modes can compete with the two-gluon channel, squarks should not be 
dramatically heavier than gluinos. On the other hand, to avoid dominance of the single- 
gluino decay, the decay mode g —>■ qq should be kinematically forbidden, corresponding to 
m g < m q- We therefore evaluate the corresponding decay rates and study if, given suitable 
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Table 1: Color interaction of two SU(3) octets. 



choices for the respective masses, the two-photon mode might lead to a possible signal. 
In Section [6] a rough estimate of the dominant signal versus background is presented. 
Section [7] contains our conclusions. 

2 Properties of gluinonia 

Let us start with a brief recapitulation of the quantum numbers of gluino boundstates 
and the corresponding color, spin and orbital momentum configurations [1-3]. These 
differ from those of quarkonia due to the restrictions arising from the Majorana nature of 
gluinos, and due to their different color assignment. 

Two color-octet states can be combined into irreducible representations as follows 

8® 8 = l s ©8 s ©8 a ©10 a ©T0a©27 s , (1) 

where the index indicates the (anti-)symmetry with respect to their color index. The 
interaction can be either attractive, repulsive or absent (in lowest order). In lowest order 
the coefficient of the QCD potential is given by the expectation value of the product of the 
color generators F^F^, taken between two-particle states in the respective representation. 

This product, in turn, can be expressed by the eigenvalues of the quadratic Casimir 
operator of the constituents, Ca = 3, and the boundstate in representation R, Cr = 
(F«) 2 : 

F ^ . F a > 2 = I [( F R ) 2 - (F"' 1 ) 2 - (F a - 2 ) 2 ] = 1 (C R - 2C A ) . (2) 

The results are listed in Tab. [U In the following we shall limit the discussion to the cases 
with negative coefficients, corresponding to attraction. 

Fermi statistics and the Majorana nature of the gluinos lead to additional restrictions. 
For the symmetric color configurations l s and 8 S antisymmetric spin-angular momentum 
wave functions, (— 1) L+S = 1, are required, for the antisymmetric color configuration 8 a 
symmetric ones, (— 1) L+S = —1. The intrinsic parity of a Majorana particle can be chosen 
to be imaginary, its parity under charge reflection real, leading to negative intrinsic parity 
and positive charge parity of the boundstate. For a few lowest orbital angular momenta 
the complete set of boundstate quantum numbers is listed in Tab. El 
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Table 2: Lowest-lying states J PC of the gluinonium spectrum. 



The production rate of non-relativistic boundstates in hard collisions is proportional to 
the squared wave function at the origin (for S waves) or its derivative (for P waves), 
the latter being significantly suppressed and of relative order v 2 . For this reason the 
following discussion will be limited to S waves only. We, furthermore, anticipate that the 
experimental mass resolution and, eventually, the large natural width of the boundstate 
will lead to a sizable smearing of the "narrow" resonances. The signal to background ratio 
will be small and their detection difficult. For this reason we will consider color-singlet 
states only, where the level spacings and production rates are enhanced approximately by 
the square of the color coefficient listed in Tab. [TJ 

Let us now discuss the QCD potential. Similarly to the case of heavy-quark boundstates 
it can be decomposed into the following two terms 



which are ordered according to their inverse powers in the constituent mass m and are 
given in coordinate space in close similarity to the quarkonium potential (see e.g. [19]). 
For S waves it reads 
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Here a s stands for the strong coupling constant in the MS scheme, nj is the number of 
active quark flavors, Ca = 3, C F = 4/3, Tp = 1/2, 7^ = 0.5772 . . . the Euler constant, 
C(3) = 1.2021 . . . and rrig is the constituent pole mass. 

For the Coulomb part Vc the transformation from color-triplet, where the result can 
be found in the literature, to constituents in an arbitrary representation R, bound to a 
singlet state, can be understood as follows: In LO the overall normalization is changed 
by substituting Cp — » Cr. This is the only modification also in NLO, thus leaving the 
constant a\ and the terms proportional to (3q unchanged. For the case under discussion 
(R = A) this is evident from Fig. [TJ Diagram (a) contributes proportional C R , diagram 
(b) proportional C\ — \CrCa- Its C R part can be combined with (a) and in total the 
ladder plus crossed ladder can be collected in the Cr/t part of the potential. The CrCa 
part of diagram (b) contributes to the Ca term in f3 and in the constant a\. Diagram 
(d) contributes proportional CpTprij and is responsible for the T F Uf terms in (3 and a±. 
The corresponding considerations are also applicable to the two-loop case, see e.g. Fig. [21 
For diagrams with virtual gluons only (e.g. Fig. [2] (a)), the substitution Cp —> Cr is 
valid throughout. For diagrams with exactly one dressed gluon propagator exchange (e.g. 
Fig. [2] (c), (d), (e)), the same substitution is valid, as far as the coupling of the gluon to 
the constituent is concerned and only the external Cp is replaced by Cr. Diagrams with 
three-gluon coupling through a triangular fermion loop (Fig. [2] (b)) are also proportional 
to CrCaTfUj and it is again the first factor, which is chosen to be C F in the case of quark- 
and Ca in the case of gluino-constituents. Diagrams similar to those of Fig. CE] (a)-(c), but 
with one dressed gluon (e.g. Fig. [2] (f)) can be handled like the NLO diagrams. In total, 
up to NNLO, the Coulomb part of the potential for the binding of constituents in an 
arbitrary representation R is obtained through the substitution Cp — * Cr in the overall 
fact or 

The term U3 of Eq. (T4J) of order a 2 s /m originates from a one- loop calculation, and the 
factor C F (C F /2 — Ca) must be replaced by Cr(Cr/2 — Ca)- The terms of order a s /m 2 



1 Wc do not expect that this statement holds true to arbitrary high orders. In contrast to heavy quarks 
gluinos together with gluons can form color-singlet states and even in quenched QCD the static potential 
between two gluinos is not expected to be confining. 
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Figure 1: Typical NLO contributions to the qq potential. Straight and curely lines repre- 
sent quarks and gluons, respectively 




Figure 2: Typical NNLO contributions to the qq potential. The same notation as in Fig. [I] 
has been adopted. 



collected in U\ and U 2 can be translated simply by replacing the global factor Cp by Cr. 
For our application the substitution Cr —> Ca is implicit. 



3 Spectroscopy 

Using the potential as specified in Eqs. and (J4]), binding energies, level spacings and 
bound-state wave functions are easily obtained by solving Schrodinger's equation numer- 
ically for Vc and adding the singular terms collected in Ui, U<i and U3 in perturbation 
theory or, alternatively, by evaluating energy level and wave functions in a perturbative 
series in a s . 

The analytic results both for binding energies E n and wave functions at the origin \^f n (0) | 2 
are listed in Appendix In numerical form the results for the two lowest levels are given 
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by 




{l + a s (2.44Li + 3.20) 



+a 2 s [(4.47L? + 9.71Li + 12.47) c + (20.81) nC ] } , 

E 2 = -H^M { i + as (2 .44L 2 + 4.42) 
lb 

W s [{4A7L 2 2 + 14.19L 2 + 20.54) c + (11.95) nC ] } , 

*i(0)| = g Q {1 + ^(3.66^-0.43) 

+a* [(8.93L? - 5.11Li + 5.83) c + (27Li + 57.38) nC7 ] } , 

^2 (0)| 2 = ° A {1 + a, (3.66L 2 - 0.18) 



with L n = In [n\ij [rrigC a^s)) and being the gluino pole mass. The first of the a 2 terms 
give the Coulombic corrections (C) the second ones the non-Coulombic ones (nC). Here 
and below we use a s = a s (/x) as defined in the MS scheme and the /1 dependence from 
the numerical solution of the three loop renormalization group equation with the starting 
value a s {Mz) = 0.1176. We use n/ = 5. This is well justified, as long as the characteristic 
scale a s rrig is smaller than the top-quark mass, an assumption well justified for all gluino 
masses under consideration. Top-quark effects could be included following [37,38]. If 
not stated otherwise, in the following we use for the renormalization scale \i — fig with 
Us = rrigCAOt s {iis)/ n corresponding to L n = 0. 

We have convinced ourselves that the perturbative series and the numerical results are in 
very good agreement as far as Vc is concerned. The subsequent analysis will therefore be 
based on the fully perturbative approach. 

The groundstate energy E\ is shown in Fig. [3] as a function of 2m,g for gluino masses rrig 
between 300 GeV and 1500 GeV. The LO (dotted), NLO (dashed) and NNLO predictions 
are displayed individually. For the NNLO result both the result for Vc (dash-dotted) 
and the one including the corrections collected in Ui, U2 and U3 (solid) are shown. As 
stated above, nig is understood as the constituent pole mass. The poor convergence of 
the perturbative series can partly be traced back to the large non-Coulombic correction 
which appear for the first time in NNLO and amount to twice the NNLO Coulombic one. 
Qualitatively this behavior is quite similar to the one observed for quarkonia in [25,39,40]. 
It remains to be seen to which extent inclusion of orders of N 3 LO and higher will stabilize 
these predictions]! In Fig. H] we show as an alternative representation the predictions for 
the ground state energy, using as reference the potential subtracted mass [22], which is 

2 For the quarkonium case, see Refs. [25,41]. 
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Figure 3: Ground state energy E\ in the pole mass scheme as function 
of twice the constituent mass. 



related to the pole mass through 
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In the following we adopt for the factorization scale [if = rrigCACtsins), independent of 
n (with us = m g CAOi s (ns) l n as before). (Note, that mps now also depends on n.) We 
refrain from listing explicitly the results for the energy levels in the potential subtracted 
scheme since the corresponding formulae are quite bulky. They can easily be obtained 
from Eqs. (jSJ) and (JTj). 

The implementation of the potential subtracted mass leads to a significantly improved 
convergence of the perturbative series for Ef s as far as the Coulombic part is concerned 
(see Fig. HJ) and similarly for E^. However, the meson mass difference between the 15 
and the 25* state, M(2S) — M (IS) , is independent of the choice of the scheme and exhibits 
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Figure 4: Ground state energy E\ in the potential subtracted mass 
scheme as function of twice the potential subtracted mass. The curves 
have been obtained using the pole mass as input and evaluating both the 
potential subtracted mass (using Eq. (ED) and E[ s to a given order in 
a s . Note that the dash-dotted and dashed curves are almost on top of 
each other. 



a poor convergence in both cases. The predictions for M(2S) — M(1S) are shown in Fig. [5l 
adopting fi = fi$- The NNLO terms are evidently important and the poor convergence of 
the prediction for E\ is reflected in this figure. 

In the on-shell scheme the 15* binding energy is evidently closely related to the energy 
difference between the groundstate and the onset of open gluino production or, more 
specifically, the threshold for pair production of color-neutral (gg) hadrons. However, 
anticipating a mass resolution of several tens of GeV at least, it will be difficult to resolve 
the densely distributed radial excitations. For this reason the energy difference between 
the groundstate and the first radial excitation is a convenient measure of the isolation of 
the IS state. 

Production and decay of non-relativistic S-w&ve bound states are proportional to the 
square of the wave function at the origin. The analytic results are collected in the Ap- 
pendix |AJ numerical results are listed in Eq. and shown in Fig.EJ The enormous size in 
particular of the non-Coulombic NNLO terms has also been observed for the top-antitop 
system and destabilizes the predictions [32]. At present the size of these corrections must 
be considered as an estimate of the theory uncertainty. 
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Figure 5: Meson mass difference M(2S) — M(1S) as function of twice 
the constituent mass. 



The renormalization scale dependence of the LO, NLO and NNLOc results for E\ and 
|*i (0)r are shown in Figs. [7] and [8] for a constituent mass of 1 TeV. The predictions are 
normalized to the LO prediction evaluated at the scale /x as defined above. A significant 
stabilization is observed with increasing order in a s and, furthermore, the higher order 
contributions are reasonably small for the preferred choice of the renormalization scale, 
at least, as far as the Coulombic part is concerned. This is in particular true for the 
ground state energy parameterized in terms of the potential subtracted mass where close 
to \x = 150 GeV the Coulombic NNLO results even vanish. In the case of the wave 
function one observes around [i = 300 GeV corrections of the order of a few percent from 
NLO and NNLOc, however, a huge contribution from the non-Coulombic terms. 

4 Bound state production and decay 

In the present context squarks are assumed to be significantly heavier than gluinos, such 
that the direct decay g — > qg is forbidden. An extreme example is provided by Split SUSY 
[7, 8] with nig ^> m~ g and we shall adopt this simplifying mass assignment throughout this 
section. However, this assumption could be significantly relaxed. Indeed, once rrig > rrig 
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Figure 6: Squared ground state wave function at the origin as function 
of twice the constituent mass. 



is valid the gluinos have a very small decay rate, for rriq 3> m g of order [42] 
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characterizing the rate for gluino decays into neutralinos (similar estimates being applica- 
ble for decays into charginos) such that the decay rate is far smaller than the annihilation 

decay 

The evaluation of the rate for gluinonium annihilation into gluon (and quark) jets proceeds 
similar to the one for quarkonia. In Born approximation the color factor is replaced as 
follows 
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and another relative factor \ arises from the Majorana nature of the gluinos, whence [1-3] 
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with |i?(0)| 2 = 47r|\l/ 1 (0)| 2 . Using the wave function obtained in lowest order perturbation 
theory leads to T gg rs {CAOi s ) b mg/A which provides a qualitative estimate of the full result. 
The perturbative corrections arising from virtual and real emission can be calculated 
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similar to the ones for quarkonia [43]. They reduce the scale dependence of the lowest 
order approximation and are complementary to the higher order corrections to the wave 
function. In total one finds 



NLO 
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n f T F 



(11) 

Here a s is understood in the theory with rif = 6 effective flavors and all quark masses, 
including m t , have been neglected. The term present in |i?(0)| 2 is interpreted in the 
rif — 5 effective theory. The predictions for IYo and Tnlo are shown in Fig. [91 where the 
NNLO wave function as shown in Fig. [6] (excluding the non-Coulombic contribution) and 
fi = 2rrig for the renormalization scale in Eq. ( TTTj) have been adopted. As anticipated, the 
decay rate of the bound state is small compared to the level spacing and large compared to 
T(g — > ^qq), the decay rate of a single gluino, once rrig < rrig. Hence annihilation decays 
will constitute the dominant signal. However, it should be emphasized, that bound states 
will exist even for m g > rrig, as long as T (g — > qq) < E2 — E\, which is true for the region 
nig ^ 0.9 ms. 
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Figure 9: Decay rate of the IS state into gluons. 



The bound state production cross section can be calculated similarly to the one for the 
pseudoscalar bound state of two top quarks (usually called r) t ) [43]. In Born approximation 
the reaction proceeds through gluon fusion and the radiative corrections involve additional 
gq, gq and qq initiated subprocesses. The final result is quite similar to the one for 
quarkonia production [43]. We find 
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Here z 

again defined in the nf = 6 flavor theory, are the Altarelli-Parisi splitting functions and 
[...]+ denotes the usual + distribution. For the wave function at the origin we again adopt 
the NNLO result without the non-Coulombic terms. The predictions for the production 
cross section are shown in Fig. [TU1 where the parton distribution functions MSTW2008LO 
(MSTW2008NLO) [44] have been used for the LO (NLO) calculation. Leading order 
(dashed) and NLO (solid curve) predictions are evaluated using /i = \i F = 2rrig. 

For a gluino mass of 1 TeV the /i dependence (with \i F — /i) of the LO and the NLO 
results is shown in Fig. [TTJ together with the decomposition of the NLO result into the 
contributions from the different subprocesses. 

From Figs. [10] and [11] we learn that the NLO correction amounts to more than +70% 
fairly independent of the renormalization scale. The dominant uncertainty, not displayed 
in these figures, arises from the wave function at the origin, with indications for an increase 
by another factor of two. 

For an anticipated integrated luminosity of 100 fb _1 between 10 and 10 6 events are ex- 
pected for rrig between 1.5 TeV and 300 GeV. It remains to be seen, if these can be 
separated from the hadronic continuum background. 



5 Annihilation decays into top quarks and into two 
photons 

In the limit of extremely heavy squarks (and conserved R-parity), corresponding to Split 
SUSY, decay modes of single gluinos are extremely suppressed, such that gluinos (con- 
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Figure 10: LO and NLO production cross section for the IS state at 
y/s = 14 TeV. 



fined in color-neutral hadrons like gg) could even travel microscopic distances. In this 
case gluinonium annihilation into two gluons is the only relevant channel. However, once 
squark and gluino masses are comparable, new annihilation decays become possible. At 
first glance one might consider the mode (gg) —>■ qq, mediated by squark exchange. How- 
ever, for massless quarks this process vanishes as a consequence of helicity conservation. 
For massive top quarks one finds 



using mi = = m^. The complete result for m~ tx ^ and including squark mixing 
is given in Appendix IB. II 

The result is plotted in Fig. [T2] for different ratios of m^/rrig. Branching ratios of up to 
7% are possible for gluino masses close to 300 GeV and squark masses not much bigger. 
For heavier squark and gluino masses the ratio drops quickly. 

As a second possibility we consider the loop induced decay into two photons. Although 
suppressed by the electromagnetic couplings a 2 , the cleaner signal-to-background ratio 
might (at least in principle) compensate this disadvantage. For arbitrary squark mixing 
and masses the result is given in Appendix IB.2I In the limit of massless quarks and 
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Figure 11: Renormalization- and factorization-scale dependence of the 
IS* production cross section and decomposition into contributions from 
various subprocesses. The line labeled "gg" contains also the LO result. 



degenerate squarks the decay rate simplifies considerably 

2 
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With a ratio i? 77 around one to several 10 5 the two-photon signal will not be detectable. 



6 Signal versus background 

Assuming the two-gluon mode to dominate the decay of gluino bound states, the resonant 
signal must be isolated from the two-jet continuum. Below we present an order of magni- 
tude estimate for the signal-to-background ratio, following essentially the arguments given 
in Ref. [3]. For simplicity it is assumed that quark and gluon jets can be experimentally 
separated. In this case the background is dominated by gluon-gluon scattering. Further- 
more, we enhance the signal-to-background ratio by excluding small scattering angles in 
the partonic center-of-mass system, requiring |cos(0)| < z < 1. 

To isolate a resonance in the gluon-gluon channel, a good experimental resolution AM 
of the invariant dijet mass M,-,- is essential. Let us compare the differential cross section, 
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integrated over the interval M — < Mjj < M + for background and signal, and 
integrated over the same angular range. This leads to the signal-to-background ratio 



r(z, AM) 



S(z, AM) _f r 
B(z, AM) 



3 AMa 2 s 



129 - 32z 2 - z 4 



1 - 



-.2 



24, fl + z 
— log 

z \ 1 — z 



(16) 



Using AM/Mjj = 0.038 + 38/Mjj for the dijet resolution [45], one obtains the signal-to- 
background ratio as shown in Fig. [T31 varying between 0.4% to 0.7%. Given event rates 
between 10 3 and 10 5 for rrig between 500 GeV and 1 TeV the detection of gluino bound 
states might become feasible, in particular at a high luminosity version of the LHC. 



7 Conclusions 

For a large class of supersymmetric extensions of the Standard Model gluinos will be 
sufficiently stable to form color-neutral non-relativistic boundstates. The NNLO potential 
between color-octet constituents forming color singlets has been derived. Just as in leading 
and next-to-leading order the Coulombic part of the NNLO potential is related to the 
quarkonium potential through the substitution Cp — > Ca in the overall factor. The 
transformation of the short distance terms involving an additional factor 1/m is only 
slightly more complicated and discussed in Section [2j Based on this potential, predictions 
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Figure 13: Signal-to-background ratio for different values of z. 



for the binding energy and the wave function at the origin of the lowest-lying boundstates 
and the first radial excitation have been derived. For gluino masses between 300 GeV 
and 1500 GeV the IS* binding energy (in the pole mass scheme) varies between 15 GeV 
and 50 GeV, the IS — 2S energy difference between 10 GeV and 30 GeV. For the case of 
gluinos being lighter than squarks the gluino decay rate is significantly smaller than the 
level spacing and the boundstates are thus well separated from the continuum. Also the 
wave function at the origin, which enters bound state production and decay rates, has 
been evaluated in NNLO. For m g > rrig the decay rate of gluinonia into two gluon jets is 
evaluated in NLO and shown to dominate all other modes. For top squarks with masses 
comparable to m g the decay mode into it may reach several percent, the decay through 
virtual quark loops into two photons is too small to be detectable. We evaluate the 
boundstate production cross section again in NLO. Large positive corrections of around 
50% arise from the NNLO terms in the wave function at the origin and around 70% 
from the NLO terms in the cross section. The production rate depends strongly on the 
gluino mass, with up to 10 6 events for gluino masses close to the lower limit around 
300 GeV (for an integrated luminosity of 100 fb _1 ) dropping quickly to 10 events only for 
rrig = 1.5 TeV. Using standard assumptions for the dijet-mass resolution and assuming 
that gluon and quark jets can be separated, a signal-to-background ratio between 0.4 and 
0.7% percent might be conceivable. For favorable gluino- and squark-mass combinations 
gluinonia might thus be detectable. 
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A Energy levels and wave functions 



For completeness we present the perturbative analytic results for the energy levels E n 
of the S waves and the corresponding values of the squared wave function at the origin 



', truncating the perturbative series at NNLO. We present the results following 



the notation of [39]. We define 

En = 

l^(o)l 2 = l^ 0) 

and parameterize the corrections as 
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where C stands for the corrections to the Coulombic part of the potential and nC for the 
remaining one. The LO values are given as 
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The NLO and NNLO corrections are 
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The constants ai,a 2 ,0o and 0i depend on the Casimir operators Ca,Cf and T F and are 
defined in Eq. (JSJ). The harmonic sums Si(n), the nested harmonic sums Si t j(n) and the 
zeta-function £(z) are given as 
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Figure 14: Feynman diagrams for the decay of (gg) into tt. 



B Amplitudes for the decays into top quarks and 
into photons 

B.l Decay into top quarks 

In Fig. [TH the corresponding diagrams for the decay of two gluinos into top-antitop are 
shown. Diagram (a) drops out for the case of interest because the two gluinos are projected 
here onto a color-singlet state. The resulting amplitude for (j3) and (7) reads 
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where si,S2,s and s are the spin indices of the top, the anti-top and of the gluinos, P_ 
and P + are the left- and right-handed projectors as usual and the Um n are the elements 
of the orthogonal squark-mixing matrix 
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(23) 



The relative minus sign between the diagrams {(3) and (7) of Fig. [T3] stems from the 
interchange of two fermionic operators. 

The amplitude in Eq. (1221) has now, following the methods explained in Refs. [47, 48] , to 
be projected onto the spin-0 state using the Clebsch-Gordon coefficients (0, 0; s,s), where 
the squared brackets of Eq. ( 1221) become due to the antisymmetry of the spin function 
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Figure 15: Representative diagrams for the decay into two photons at LO. 



Now the calculation can be done as for the case of the decay into two gluons and the rate 
reads 

2 m- (m+ — 2m~[/®C/® N " 



B.2 Decay into photons 

Treating all the virtual quarks in Fig. [15] as massless, the onto spin and color-projected 
amplitude has the form 




A((gg)i s -> 77) 
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(26) 



The sum in the first equation of (I26p goes over the flavors (/) and the indices of the 
squarks (h) . The ki are the momenta of the outgoing photons, the 6j their polarization 
vectors and the Qf the charges of the corresponding quarks and squarks. The definition 
of the Co -functions is the same as in Ref. [49]. 
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